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ABSTRACT 

Context. Even if the '-C+'"C reaction plays a central role in the ignition of Type la supernovae (SNIa), the experimental determination 
of its cross-section at astrophysically relevant energies {E < 2 MeV) has never been made. The profusion of resonances throughout 
the measured energy range has led to speculation that there is an unknown resonance at £0 ~ 1.5 MeV possibly as strong as the one 
measured for the resonance at 2.14 MeV, i.e. (wy)^ = 0.13 meV 

Aims. We study the implications that such a resonance would have for our knowledge of the physics of SNIa, paying special attention 
to the phases that go from the crossing of the ignition curve to the dynamical event. 

Methods. We use one-dimensional hydrostatic and hydrodynamic codes to follow the evolution of accreting white dwarfs until they 
grow close to the Chandrasekhar mass and explode as SNIa. In our simulations, we account for a low-energy resonance by exploring 
the parameter space allowed by experimental data. 

Results. A change in the '-C+'^C rate similar to the one explored here would have profound consequences for the physical conditions 
in the SNIa explosion, namely the central density, neutronization, thermal profile, mass of the convective core, location of the runaway 
hot spot, or time elapsed since crossing the ignition curve. For instance, with the largest resonance strength we use, the time elapsed 
since crossing the ignition curve to the supernova event is shorter by a factor ten than for models using the standard rate of '-C+'~C, 
and the runaway temperature is reduced from ~ 8.14 x 10^ K to ~ 4.26 x 10^ K. On the other hand, a resonance at 1.5 MeV, with 
a strength ten thousand times smaller than the one measured at 2.14 MeV, but with an a/p yield ratio substantially different from 1 
would have a sizeable impact on the degree of neutronization of matter during carbon simmering. 

Conclusions. A robust understanding of the links between SNIa properties and their progenitors will not be attained until the '^C-l-'^C 
reaction rate is measured at energies ~ 1.5 MeV. 

Key words, nuclear reactions, nucleosynthesis, abundances - supernovae: general - white dwarfs 



[ 1. Introduction 

The idea that Type la supernovae (SNIa) are triggered by ex- 
plosive carbon burnin g in degenerate r nateria l dates back to 
, the seminal paper of iHovle & Fowleil (Il960h : ' 'Pure carbon 
, is explosive through ^^CC^C,af°Ne and i2c(i2c,p)23Na at a 
temperature somewhat less than 1.5 x 10^ degrees, even on a 
timescale as short as 1 second". Since that time, the paradigm 
of SNIa has changed fro m carbon detonation (lArnettI [19691) . 
first, to carbon deflagration (iNomoto et al.lll97"6l) . and, finally, to 
the cuiTently favoured paradigm: a delayed detonation starting 
as a carbon deflagration (Khokhlov 1991, delay ed-detonation 
model). S hortly after th e form ulati on of the delayed-detonation 
paradigm. lH"oflichetalJ (Il995h and lHoflich & Khokhlov! d 19961) 
showed that this kind of model can account for the basic obser- 
vational properties (brightness, light curves) of normal and sub- 
luminous SNIa. Nowadays, sofisticated three-dimensional sim- 
ulations of SNIa base d on different flavours of the delayed- 
detonation model (e.g. iGamezo et allllool: iMeakin et alJl2009t 



IBravoet all 120091: iRopkeet all 1201 lb are confronted with data 
in order to explain i ncreasingly subtle observational detail s 
(iMazzah et alJ boOTi: iKasenetalJ 120091: iMaeda et al.1 l2010l) . 
These numerical simulations have taught us that the outcome 
of the explosion is very sensitive to the runaway conditions (e.g. 
iSeitenzahl et al.ll201 ih . which are set during the so-called car- 
bon simmering phase, the pre-explosive period that begins when 
carbon starts burning at a slow pace and ends when the nuclear 
timescale becomes comparable to the white dwarf (WD) sound 
crossing time, fsound - 0.1 s (Piro & Chang 2008). The evolu- 
tion of the progenitor WD during the simmering phase (which 
is difficult to address because it is on the borderline between 
hydrostatic and hydrodynamic phenomena) is controlled by the 
competition between cooling processes (neutrino emission and 
convection) and heating by a dominant nuclear reaction: the fu- 
sion of two '^C nuclei at temperatures in the range 2 x 10^ K- 
10^ K. 
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The experimental measurement of the cross-section of the 
12q _|_i2 q reaction (carbon fusion reaction) at low energies has 
advanced slowly in the past few decades. At the timejrf formu- 
lation of the carbon detonation model for SNIa dArnetti 119691) . 
the latest known experimental results ex tended down to acen- 
ter of mass energy E^m - 3.23 MeV jPatterson et alj [T969I) . 
Seven years later, when the deflagration model was proposed, 
the ref erence '^C -h'^ C reaction rate evaluated by Fowler et aL, 
( IT971 hereafter FCZ75) was based on experimental m easure- 
ments at E c-m > 2 .45 MeV (iMazarakis & Stephenll973h . When 
[Khokhlo^ (1199 lb formulated the delayed detonation paradigm, 
the calculation of t he car bon fusion rate was based on the 
ICau^hlan & Fowleil dHM hereafter CF88) fit to experimen- 
tal cross-sections down to the same e nergy Egg (~ 2.5 M eV), 
although with higher resolution (e.g. iBecker et al.lll98l"l) . The 
rate from CF88 differs from that of FCZ75 by less than 20% 
at r > 10'' K, and by less than a factor 2.5 at T > 10** K. 
In the past 30 years, the lower limit to E^m reached experi- 
mentally has decreased to only 2.10 MeV (lAguilera et"a 
iBarron-Palos et all l2006t ISpiUane et all |2007), because of ex- 
perimental difficulties in reducing the background related to 
secondary reactions induced by hydrogen and deuterium con- 
tamination in the carbon targets. Present-day data exhibit pro- 
nounced resonances throughout the measured energy range. 
ISpillane et al ] (I2007h . indeed, found a strong resonance at /sr - 
2.14 MeV, close to the low-energy limit of their measurements. 
It turns out that the Gamow energies relevant to the carbon fu- 
sion reaction during the simmering process lie far below the 
lower energy tested experimentally. For instance, for tempera- 
tures in the range 2x10** K- 10^ K, the Gamow energies go from 
Eg + AEc = 0.82 + 0.14 MeV to Eg + AEg = 2.4 + 0.5 MeV. 
Even though the most recent low-energy measurements give 
an average astrophysical factor ~ 2-4 times smaller than the 
value recommended by CF88, the cross-sections measured at 
Ecm ^ 3 MeV are rather uncertain, which explains why the 
i2q +12 reaction rate formulated by CF88 is still nowadays 
the reference rate for astrophysical applications. 

There has been lively discussion about the extrapolation of 
the '^C -h'^ C rate throughout the Gamow energy range. On the 
one hand, theoretical phenomenological models predicting the 
behavior of the non-resonant astrophysical factor at low ener- 
gies, based on different approaches, diff'e r by up to ~ 2-3 or- 
ders of magnitude at Ecm ~ 1-1-5 Me V (iGasques et alJl2005t 
[Aguilera et al. 2006; Jiang et al.ll2007h . Most of fliese predic- 
tions give rates well bel ow the standard rate of CF88 (see, e.g.. 
Fig. 3 in IStriede^ 12008 ). On the other hand, the profusion of 
measured resonances has led to speculation about the hypotheti- 
cal presence of a resonance within the Gamow energy range, say 
at Ecm ~ 1-5 MeV, that could completely dominate the reaction 
rate at the densities and temperatures characteristic of the carbon 
simmering phase of SNIa. 

Our aim is to explore the consequences that a resonance in 
the '^C C reaction near E^m ~ 1-5 MeV would have for 
the physics of SNIa. As a strong resonance at this energy is ex- 
pected to dominate the rate, we can disregard the uncertainties 
in the behavior of the non-resonant part of the astrophysical fac- 
tor at low energies, and adopt for it the CF88 rate. This kind 
of speculation about a hypothetical low-energy resonance in the 
carbon fusion reacti on rate has been addressed in several works. 
ICooper et al.l (l2009h studied the changes induced by such a reso- 
nance on the physics of superburst ignition on accreting neutron 
star crusts. Bennett et al. (2010a) and Bennett et al. (2010b) ex- 
plored the impact on the s-process induced by an increase in the 
carbon fusion reaction rate up to a factor 50,000 at T = 5 x 10** K. 



Finally, llapichino & Lesaffi3 (l2010l) analyzed the influence of 
uncertainties in the '^C C reaction at T < 5 x 10^ K on the 
ignition of the WD core. They found that the ignition density de- 
pends slightly on the presence of a low-energy resonance while 
the ignition temperature is almost unaffected. However, their re- 
sults are influenced by the small value of the resonance strength 
they assumed, which implied that their reaction rate exceeded 
that from CF88 by only a factor of two. 

As there exists no theoretical framework that allows us 
to predict the location and strength of the resonances in the 
i2q +12 reaction, the best approach to studying their conse- 
quences is an exploration of the parameter space of a low-energy 
resonance. In the next section, we determine the resonance en- 
ergy and strength that are compatible with the available cross- 
section data, and discuss the expected implications for our un- 
derstanding of SNIa. In Section [3] we present the results of 
our simulations of the evolution of a WD from mass accretion 
to supernova explosion. These simulations have been obtained 
with the use of one-dimensional hydrostatic, hydrodynamic, and 
nucleosynthetic codes where the '^C C reaction rate was 
computed consistently assuming different resonance properties. 
To cover a wider range of evolutionary scenarios (in particu- 
lar, a wider range of accretion rates that would lead to different 
density-temperature tracks), we adopt in Section |4] a somewhat 
different approach: we follow the evolution of a homogeneous 
region through the latest stages of accretion until the dynamical 
instability with different physico-chemical parameters. In this 
section, we study the dependence of the global neutronization of 
the convective core of the WD during carbon simmering, and the 
runaway temperature of convective bubbles. Finally, we present 
our conclusions in Section|5] 

2. Low-energy resonance 

We explore the consequences of a resonance in the neighborhood 
of £0 = 1.5 + 0.2 MeV with a strength limited by the measured 
cross-sections at low energy. The low-energy measurements that 
have provided the tightest constraints to date are those obtained 
by Spillane et al. (2007, hereafter S07), who explored energies 
£■ > 2.10 MeV and found that the resonance structure continues 
down to this energy limit. They also found a resonance at E^m = 
2.14 MeV, although this resonance has a quite limited impact on 
carbon burning: it speeds up slightly the carbon fusion rate at 
temperatures in the range T ~ 6 x 10'*-1.2 x 10^ K (e.g., by a 
factor 1.8 at r - 8 x 10** K). On the other hand, the ff-channel 
yield of this resonance is larger than the p-channel by a factor 
~ 5, at variance with the standard formulation of CF88, in which 
both channels have similar strength. 

Our '^C -h'^ C reaction rate is given by the sum of the non- 
resonant contribution (that we calculate following CF88) plus a 
resonant one that accounts for the resonance found by S07 at 
Ecm - 2.14 MeV, and the assumed low-energy ghost resonance 
(LER). We compute the contribution of both resonances to the 
carbon fusion reaction rate as (see, for instance, S07): 



{o-vf = 2.28 x 10"^'*r„ 



exp(-24.8/r9) 



"0.13 meV 



exp(-11.6£'R/r9) 



1 -1 
cm s 



(1) 



where E^ is the energy, in MeV, at which a resonance is as- 
sumed, and (ciry)j!^ is the ghost resonance strength. 

Although there exist in the literature some discussions of the 
properties of such a resonance that are compatible with current 



2 



E. Bravo et al.: SNIa and the '^C+'^C rate 



experimental data, different authors allow for different Umits. 
Hence, we repeat here the derivation of the maximum possible 
strength of a LER (e.g. Cooper et al. 2009). 

As a first approach, we demand that the ghost resonance at 
/sr contributes to the cross-section at E^m = 2.10 MeV less 
than 10% of the value measured by S07 at the same energy, 
crexp(2.1 MeV) < 0.8 nb. Assuming a narrow resonance, as are 
all known resonances of the carbon fusion reaction, Fr 40- 
100 keV dAguilera et al.l ^06), the energy de pendent cross- 
section is given by t he Breit-Wigner formula (e.g. IClavtonI 1 968t 
iRolfs & Rodnevl[l9 88): 



crR(E) 



0.6566 



((yy)R Tr 



Ae (MeV) (E - ErY + rl/4- 



(2) 



where A = 6 is the reduced baryon number Substituting E 
2.10 MeV, Tr = 100 keV, and ctr < 0.08 nb we get 



(cjj\ < 7.3 X 10""'£(MeV) 



(E-ERf+Fl/A 
Tr 



= 15.4 [(2.10 -£r)^ +2.5 X 10"^] meV, 



(3) 



where Er is in MeV. Thus, the resonance strength is limited to be 
less than 9.9, 5.6, and 2.5 meV for £r = 1.3, 1.5, and 1.7 MeV, 
respectively. We note that a smaller value of the resonance width, 
Fr < 100 keV, would result in still higher upper limits to the 
resonance strength. 

A second, similar, approach to the determination of the max- 
imum possible strength of the ghost resonance makes use of the 
recommended value of the modified astrophysical factor, S *, de- 
fined as 



S'(E)^cr (E) E exp (87.21£"'^2 + 0.46£) , 



(4) 



where E is in MeV. The averag e experirnental 5 * at the lowest 
energies is ~ 0.7 x 10^^ MeV b (Strieder'2008'), while the value 
recommended by CF88 is 3 x lO"" MeV b. Taking an approx- 
imate value of the experimental modified astrophysical factor 
5exp (2.10 MeV) < 10'^ MeV b, and using Eq.|2]for the energy- 
dependent cross-section due to the ghost resonance, we obtain 



(w7)r < 



(E - £r)2 + r2 /4 



0.6566 MeV bai-n Fr 
xexp(-87.21£"'^2-0.46£) , 

which, after substituting £ = 2.10 MeV and Fr - 100 keV, gives 



(5) 



(wy)R < 2.54 [(2.10 - Er)^ + 2.5 x 10"^] meV , 



(6) 



again with Er in MeV. In this case, the resonance strength is lim- 
ited to 1.6, 0.92, and 0.41 meV for £r = 1.3, 1.5, and 1.7 MeV, 
respectively. These limits increase by a factor ~ 10 if the reso- 
nance width is Fr = 10 keV. 

A quite different approach, found in the literature, consists 
in assuming that the modified astrophysical factor remains con- 
stant at energies below 2 MeV. This assumption leads to much 
lower resonance strengths at the resonance energies we explore. 
However, we note that these are not upper limits to the resonance 
strength compatible with experimental data. 

In Fig. [U (top panel), we show the modified astrophysical 
factor of the '^C+'^C reaction as a function of energy assum- 
ing a strong LER, compatible with the limits derived above, 
5 * (2.10 MeV) < lO"' MeV b. As might be expected, the farther 




1.6 1. 
E (MeV) 



2.2 




0.5 1 1.5 2 

Temperature (10" K) 



2.5 



Fig. 1. Astrophysical factor and rate of the '^C+'^C reaction 
due to a resonance at an energy Er - 1.3 MeV (dashed line), 
1.5 MeV (solid line), and 1.7 MeV (dot-dashed line), with reso- 
nance strengths of 16 meV, 9.2 meV, and 4.1 meV, respectively 
Top. Modified astrophysical factor as a function of energy. The 
bar at the bottom right shows the maximal astrophysical factor 
measured by S07 at an energy of ~ 2.1 MeV (cr < 0.8 nb), which 
is the most constraining low energy measurement to date. The 
quoted strengths of the resonances would be compatible with a 
cr two orders of magnitude smaller than that measured by S07 
(if the resonance width is Fr x 10 keV). Bottom. Ratio of our 
'^C+'^C reaction, accounting for a LER, to the CF88 reaction 
rate, as a function of temperature. The dotted line close to the 
horizontal axis is the ratio of the rate accounting for the reso- 
nance found by S07 to the CF88 rate. The arrows mark the tem- 
perature ranges relevant to carbon ignition, carbon simmering, 
and the explosive burning in a subsonic flame or a detonation. 
The gap between carbon simmering and flames reflects the dom- 
inance of thermal conduction and shock heating over nuclear en- 
ergy release in stationary explosive burning fronts. 



the energy of the assumed resonance is from the measured ener- 
gies, the larger is the maximum resonance strength that is com- 
patible with the data. We can see that a shift of -0.2 MeV in Er 
allows us to increase the peak of the astrophysical factor by ap- 
proximately two orders of magnitude. The bottom panel of Fig. [T] 
shows the ratio of our '^C+'^C rate, which we calculated as the 
sum of the CF88 rate and the resonant rate given by Eq. [1] to 
the non-resonant rate (CF88), as a function of temperature. The 
assumption of a strong ghost resonance of the aforementioned 
characteristics leads to a substantial enhancement of the carbon 
fusion rate in the temperature range 2 x 10*^- 2 x 10^ K. While 
this temperature range depends scarcely on the assumed Er (the 
peak of the rate ratio shifts from 3.9 x 10'^ K to 5.8 x 10'^ K when 
the resonance energy changes from Er- 1.3 MeV to 1.7 MeV), 
the maximum ratio changes by as much as four orders of magni- 
tude between the extreme explored values of Er. 
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We can gain some insight into the impact of a ghost reso- 
nance on the physics of SNIa by analyzing the different phases 
involved in the process of explosive carbon burning: 

- At first, carbon competes with neutrino emission, until the 
temperature is high enough to ensure that nuclear heating 
cannot be balanced by neutrino cooling (ignition condition). 
In this work, we define the ignition curve as the local p-T 
conditions for which the nuclear timescale equals the neu- 
trino cooling timescale (see Fig. |2] for the definition of the 
nuclear timescale). Typically, this occurs at T ~ 2 x 10** K, 
just at the border of the temperature range where the reso- 
nance begins to noticeably increase the carbon fusion rate. 
The evolution of the white dwarf during this phase deter- 
mines the place where carbon ignites: either at the center or 
close to the white dwarf surface, the latter giving rise to an 
accretion-induced collapse of the white dwarf rather than a 
supernova event. 

- When the neutrino cooling is negligible, relative to the 
nuclear energy release, convection becomes the dominant 
energy-transport mechanism. An adiabatic thermal profile 
is imprinted in the white dwarf as long as the convective 
timescale ~ 10 - 100 s remains below the nuclear timescale. 
If the carbon fusion rate is just given by the non-resonant 
contribution, the maximum temperature at which convection 
is able to remove the excess nuclear heat is ~ 7 - 9 x 10^ K. 
When neither neutrinos nor convection can equilibrate the 
nuclear heating, a flame starts to propagate through the white 
dwarf. Here, we define the dynamical curves as the lo- 
cal p-T conditions for which the nuclear timescale equals 
the (global) convective timescale. In general, the convective 
timescale is difficult to determine and may change over time. 
However, for the purposes of studying the effects of a LER, 
it is enough to use an approximate, constant value, which we 
have taken to be t^om = 10 sQ 

Once a combustion front is present in the white dwarf, it prop- 
agates in a way that is nearly independent of the details of the 
nuclear reactions, as long as they are brought to completion. The 
front can propagate either as a subsonic flame or as a supersonic 
detonation. 

- The (microscopic) thermal structure of a subsonic flame (de- 
flagration) is determined by heat diffusion from hot ashes to 
cool fuel until the temperature of the last one is as high as 
2 - 5 X 10^ K (depending on fuel density). At higher temper- 
atures, the flux of heat released by nuclear reactions is higher 
than that of the heat diffusion. 

- Detonations are supersonic combustion waves where igni- 
tion is triggered by an initial temperature jump due to a shock 
front. In the conditions of thermonuclear supernovae, this 
jump raises the fuel temperature well above ~ 2 x 10^ K. 
Later on, the heat released by thermonuclear reactions deter- 
mines the thermal profile. 

As can be seen in Fig. [1] the range of temperatures for which 
nuclear reactions control the thermal evolution of thermonuclear 
combustion fronts in SNIa is above the range of temperatures 
at which a ghost resonance significantly affects the rate of the 

' In the literature, the term ignition is often used to describe the con- 
ditions (density, temperature, chemical composition, and geometry of 
the runaway kernels) when the white dwarf evolution becomes dynam- 
ical. This should not be confused with our use of the term ignition, re- 
ferred to in the ignition curve, which takes place tens, or even hundreds, 
of years before the actual explosion. 



Table 1. Summary of C-h C rates explored in this paper. 



Model 




i<^7k 


a/p 




(MeV) 


(meV) 




CF88 








LER-1.5-0.1-1.0 


1.5 


0.1 


1.0 


LER-1.5-0.001-1.0 


1.5 


0.001 


1.0 


LER-1.5-10.-1.0 


1.5 


10. 


1.0 


LER-1.3-0.1-1.0 


1.3 


0.1 


1.0 


LER-1.7-0.1-1.0 


1.7 


0.1 


1.0 


LER-1.5-0.1-5.0 


1.5 


0.1 


5.0 


LER-1.5-0.1-0.2 


1.5 


0.1 


0.2 



Notes. Parameters of the low-energy resonance: resonance energy, iiR, 
resonance strength, (ci>j\, and ratio of a to proton yields. 



carbon fusion reaction. Thus, we expect that the largest effect of 
such a resonance will take place during the carbon simmering 
phase. 

In this work, we systematically explore the effects of a LER 
at different values of the resonance energy, Er, and with different 
resonance strengths, (wy)^. Moreover, motivated by the exper- 
imental finding of an alpha-to-proton yield ratio, a/p, substan- 
tially different from one at the lowest measured energy (S07), we 
further study the effects of varying this ratio. As we demonstrate 
in Section l4n varying a/p has an interesting effect on the level 
of neutronization attained during the carbon simmering phase. 
In Table [T] we summarize the combinations of the resonance pa- 
rameters we have considered. Model CF88 is a reference model 
in which the '^C-h'^C rate is given just by the non-resonant CF88 
rate. The rate of all the other models is given by the sum of the 
CF88 rate and the experimentally known resonance at 2. 14 MeV 
plus a ghost resonance (Eq.[T]l. 

A first quantitative evaluation of the impact of the ghost res- 
onance parameters assumed in Table [1] can be made by com- 
paring the nuclear timescale to relevant timescales of the white 
dwarf. Figure |2] shows the nuclear timescale for the different 
models, as functions of temperature, compared to typical val- 
ues of the convective turnover timescale (fconv ~ 10-100 s). We 
see that the experimental resonance at 2.14 MeV with no ghost 
resonance (S07) has a small effect on the nuclear timescale. On 
the other hand, a resonance at Er = 1.5 MeV with a strength 
100 times smaller than that measured at 2.14 MeV (i.e., with 
(cjy)^ - 0.001 meV) can reduce the runaway temperature by 
~ 2 X 10** K, if the density is kept fixed at p = 2 X 10'' g cm"^. 
Increasing the resonance strength by two orders of magnitude 
translates into an approximate reduction in the runaway temper- 
ature by ~ 10^ K. The same reduction in the runaway temper- 
ature can be achieved by decreasing the resonance energy by 
~ 0.2 MeV. 

Similar conclusions would result from using the thermal 
timescale instead of the nuclear timescale. In general, the ther- 
mal timescale is shorter than the nuclear timescale by a fac- 
tor Tth/Tnuc - '^CpT I QNfJf\2, where Cp is the specific heat at 
constant pressure and Q is the e ffective heat release pe r carbon 
fusion reaction, Q ^ 9.1 MeV dChamulak et all 1200 8*). For in- 
stance, at r = 6 X 10^ K, Tth/Tnuc - 0.06 but, owing to the 
steep slope of the t^^^ vs. T curves (see Fig.|2|i, such a reduction 
in the timescale only slightly affects the values of the runaway 
temperatures. 

The utility of Fig. |2]is limited by the use of a fixed density, 
independent of temperature and of the parameters adopted for 
the ghost resonance. In reality, both the temperature and density 
evolve during the simmering phase and depend on the assumed 
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T (10" K) 

Fig. 2. Nuclear timescale vs. temperature. The nuclear timescale 
has been computed as t„uc = {'2-pN^Y\2 < crv where A^a is 
Avogadro's number, Yn - 0.5/12 mol g"' is the molar fraction 
of '^C, and < crv > is the rate of the '^C+'^C reaction. There are 
shown curves belonging to the CF88 non-resonant rate, the non- 
resonant plus the resonance at 2.1 MeV measured by S07, and 
the non-resonant plus the 2.1 MeV resonance plus a LER. For 
the last, the curves are labeled with the values of the resonance 
energy, in MeV, and the resonance strength, (u)y)^ in meV. 
The hatched region shows the range of the convective timescales. 
In this plot, the density is p = 2x 10^ g cm"-' and the carbon mass 
fraction is 0.5. Note that a lower carbon abundance would imply 
a longer nuclear timescale. 



'^C-h'^C rate. In the calculations reported in the next section, we 
address both dependences. 



3. White dwarf evolution from accretion to 
explosion 

We now discuss the evolution of a non-rotating WD from the 
beginning of accretion to explosion. We have followed the WD 
with an hydrostatic code until the first crossing of the dynamical 
curve at any zone within the WD. The structure so obtained was 
then fed into a hydrodynamic supernova code, where the struc- 
ture evolved until the star was completely disintegrated (max- 
imum density < 0.2 g cm"^). Both codes are one-dimensional 
and assume spherical symmetry. Strictly speaking, the onset of 
the dynamical event takes place when the evolutionary timescale 
becomes shorter than the sound crossing timescale, fsound -0.1s, 
which occurs slightly later than the crossing of the dynamical 
curve. However, when the nuclear timescale is shorter than the 
convective timescale the runaway becomes local, hence the heat 
released before attaining fsound remains small relative to the WD 
binding energy, and the WD structure does not change apprecia- 
bly in the interim. 



3.1. Methods 

Our initial model is the same as in iPiersanti et al.l (|2003h . i.e. a 
CO WD with M = 0.8 Mq, on which we accrete directly CO- 
rich matter at a rate of M = 5 x 10"^ Mq yr"'. The hydrostatic 
evolution is performed us ing the FRANEC evolutionary code 
(IChieffi & Stranierdfllll) 

with the same input physics as in 
Piersanti et al. (20031). The main difference of our approach con- 
cerns the neutrino energy losses, which are computed according 
tol^posito et al. (2002) and Esposito et al. (2003). Moreover, to 
properly describe the late part of the hydrostatic evolution, com- 
position changes produced by convective mixing are modeled by 
adding diffusion terms to the nuclear burning terms in the com- 
position equations. 

The supernova hydrodynamic c ode us ed in the present work 
is the same as in iBadenes et al.l (l2003h. in addition to both 
iBravo et al.) d 19961) and lBravoetal.1 (1201 iT) . The present mod- 
els ar e based on the delayed-detonation paradigm (Khokhlq^ 
Il991h . with a fixed deflagration-detonation transition density, 
Pddt = 1-6 X lO^gcm^^ 

3.2. Results 

The results of the hydrostatic simulation of the accretion phase 
up to the crossing of the dynamical curve are shown in Figs. [3] 
and |4] and Table [2] Fig. |3] shows the ignition and dynamical 
curves obtained using the reference CF88 rate for the '^Ch-'^C 
reaction and the changes introduced by the different sets of res- 
onance parameters that we have explored. The different ignition 
curves converge towards the CF88 ignition curve at high densi- 
ties, whereas the dynamical curves follow the opposite trend, 
i.e. they diverge from the CF88 dynamical curve as the den- 
sity increases. In all cases, as expected, the temperature defined 
by these curves, at any given density, decreases when a reso- 
nance is included. The effect, however, depends significantly on 
the resonance strength and energy. For instance, a resonance 
strength of (ajy)^ - 0.001 meV at = 1.5 MeV (model 
LER-1. 5-0. 001-1.0) scarcely affects either the ignition or the 
dynamical curves. On the other extreme, the combinations of 
((Ajy)^ - 0.1 meV at Er- 1.3 MeV, as wefl as (wy)^ = 10. meV 
at £r = 1.5 MeV (models LER- 1.3-0. 1-1.0 and LER-1.5-10.- 
1.0, respectively), result in a decrease in the temperature of the 
ignition curve by ^ 2.5 x 10** K at low densities, and a decrease 
in the temperature of the dynamical curve by ^ 4 x 10** K at high 
densities. 

The time elapsed since the beginning of accretion until the 
supernova explosion is about the same for all the '^Ch-'^C rates, 
face - 1.11-1.15 Myr, because all the models need to arrive 
close to the Chandrasekhar mass before destabilization. Models 
LER- 1.3-0. 1-1.0 and LER- 1.5- 10.- 1.0 are the ones that diff^er 
most from the reference model, for which face is approximately 
40, 000 yr shorter than that of CF88. On the other hand, the time 
elapsed between the crossing of the ignition and the dynamical 
curves (simmering time, last column of Table |2]i is quite dif- 
ferent from model to model. Even for model LER- 1.5-0. 1-1.0, 
whose strength is similar to that of the experimental resonance 
at 2.14 MeV, the simmering time is just ~ 13% of that of model 
CF88. 

Figure [3] shows the tracks followed by the runaway zone in 
the p-T plane. For models igniting at the center, the track is 
exactly the same as in CF88 up to the crossing of the ignition 
curve, as expected from the discussion in Section |2] In mod- 
els LER-1.5-0.001-1.0 and LER- 1.7-0. 1-1.0, the p-T tracks are 
hardly distinguishable from CF88, even up to the crossing of the 
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Table 2. White dwarf properties at the crossing of the ignition and dynamical curves. 



ignition curve 



Model 


r," 




X,('2c)° 






(10** K) 


(g cm-') 




(Mo) 


CF88 


2.658 


2.85 X 10'' 


0.2387 


0.5726 


LER-1. 5-0.1-1.0 


2.256 


2.56 X 10' 


0.2387 


0.5708 


LER-1. 5-0.001-1.0 


2.357 


2.84 X 10' 


0.2387 


0.5724 


LER-1. 5-10.-1.0 


2.150 


2.06 X 10' 


0.2387 


0.5678 


LER-1. 3-0.1-1.0 


1.984 


1.67 X 10' 


0.2388 


0.5627 


LER-1.7-0.1-1.0 


2.361 


2.85 X 10' 


0.2387 


0.5726 



dynamical curve 



Model 




Pc* 






M 


A/f ^ p f 


t- f 

' simm 




(10** K) 


(g cm 3) 




(Mo) 


(Mo) 


(Mo) (km) 


(yr) 


CF88 


9.715 


2.09 X 10' 


0.2237 


0.0124 


1.268 


center 


534.4 


LER-1. 5-0.1-1.0 


6.075 


2.33 X 10' 


0.2332 


0.0030 


0.996 


center 


71.6 


LER-1. 5-0.001-1.0 


7.878 


2.38 X 10' 


0.2272 


0.0078 


1.162 


center 


157.9 


LER-1. 5-10.-1.0 


5.019« 


1.90 X 10'« 


0.2708« 


0.0021 


0.889 


0.0137 149 


54.2 


LER-1. 3-0. 1-1.0 


5.240'' 


1.34 X 10"' 


0.3007'' 


0.0029 


0.996 


0.0382 233 


61.9 


LER-1.7-0.1-1.0 


7.009 


2.62 X 10' 


0.2414 


0.0050 


1.102 


center 


306.2 



Notes. Initial white dwarf mass: 0.8 Mg, accretion rate: 5 x 10"'' Mo/yr. Central temperature, density, and mass fraction and total mass of '^C 
at the crossing of the ignition curve. Central temperature, density, and mass fraction of '^C at the crossing of the dynamical curve. Total mass 
of '^C burned up to the crossing of the dynamical curve. '''' Size of the convective core. Location (mass and radius) of the first zone to cross the 
dynamical curve. Time elapsed between the crossing of the ignition and dynamical curves, or simmering time. The temperature, density and 
'-C mass fraction shown belong to the first zone to cross the dynamical curve. The corresponding values of density and '^C mass fraction at the 
center were: 2.06 x 10' g cm"' and 0.2387, respectively. '''' The temperature, density and '^C mass fraction shown belong to the first zone to cross 
the dynamical curve. The corresponding values of density and '^C mass fraction at the center were: 1.57 x 10' g cm"' and 0.2388, respectively. 



dynamical curve. However, as the dynamical curve of the ghost 
resonance models lies systematically below that of CF88, the dy- 
namical curve is reached at a different temperature and density 
than the reference model. On the other hand, model LER-1.5- 
0.1-1.0 diverges away from that of CF88 after crossing the igni- 
tion curve, which results again in both a substantially lower tem- 
perature and a slightly higher density at runaway (see Table|2|. 

The most interesting cases are those of models LER-1. 5-10.- 
1.0 and LER-1. 3-0. 1-1.0, which are characterized by either a 
large resonance strength or a low resonance energy. As can be 
seen in Tabled in both cases the dynamical curve is first reached 
off-center, at radii 149 and 233 km, respectively. The off-center 
ignition is caused by the compressional heating timescale of the 
WD being slightly longer than the thermal diffusion timescale, 
so that the WD is not completely isothermal at carbon ignition. 
A somewhat different accretion rate would provide off-center ig- 
nitions for different combinations of the resonance parameters. 
The igniting zones are located in the p-T plane to the left of the 
track of the center of the WD, and reach the dynamical curve at 
a much lower temperature and density than model CF88. Thus, 
although the presence of a LER always results in a smaller tem- 
perature at the crossing of the dynamical curve, the density can 
be either lower or higher than in model CF88 depending on the 
resonance properties. This lack of monotonicity reflects as well 
in other details of the outcome of the supernova explosion, as we 
later discuss. 

Table|2]provides additional properties of the WD at runaway, 
that are worth discussing. The size of the convective core de- 
termines the degree of dissemination of the products of carbon 
simmering throughout the WD prior to the SNIa explosion, with 
consequences for the chemical stratification of superno va rem- 
nants and the interpretation of their X-ray spectra (Bade nes et alJ 
|2008|) . In general, the mass of the convective core, Mconv, de- 
creases monotonically with increasing resonance strength or de- 



creasing resonance energy. For the extreme case, model LER- 
1.5-10.-1.0, Mconv is 30% lower than for model CF88. 

Another interesting quantity provided in Table [2] is the to- 
tal mass of ^^C that is burned during the simmering phase, 
AM('^C). This mass varies strongly from model to model, 
again monotonically decreasing with either increasing resonance 
strength or decreasing resonance energy, because of their lower 
temperature at runaway, which implies a smaller thermal con- 
tent of the WD. The value of AMC'^C) has a significant effect on 
the amount of neutronization during carbon simmering, as we 
discuss in Section l4n 

The thermal and '^C mass fraction profiles through the WD 
at the dynamical curve are shown in Fig. |4] where the range of 
the convective core can be clearly appreciated. In all the models, 
there is a temperature spike close to the surface that is caused 
by the high accretion rate, but this spike has no practical con- 
sequences for the supernova explosion. The temperature peaks 
at the center in all the models apart from LER- 1.3-0. 1-1.0 and 
LER-1. 5-10. -1.0. In these models, even though convection ex- 
tends all the way down to the very center of the WD, the '^C-h'^C 
rate increases so steeply with temperature that time-dependent, 
convective mixing is at late times only efficient in the neighbor- 
hood of the runaway zone. 

The chemical profile of '^C at the crossing of the ignition 
curve, shown as a thin line in the bottom panel of Fig. |4l re- 
flects the previous history of the WD. Above 0.8 M©, the carbon 
mass fraction is 0.5, which is the value we adopted for the ac- 
creting carbon mass fraction, while the central ~ 0.4 M© have 
a homogeneous composition with a smaller carbon mass frac- 
tion, Xc(^^C) = 0.24, as results from hydrostatic He-burning. 
At runaway, the chemical profile is the result of the competi- 
tion between convective mixing and carbon burning. When the 
dynamical curve is approached, burning consumes carbon more 
efficiently than convection can supply fresh fuel, and there is a 
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IQ7 IQR lO'" 108 10'° 

Density (g/ cm'^) 



Fig. 3. Changes in the physics of the accretion phase owing to 
a LER: ignition curves (dashed lines), dynamical curve (dot- 
dashed lines), and track of central p vs. T during the hydrostatic 
evolution (solid lines). The top left panel shows the results ob- 
tained when the '^C-h'^C rate was computed following CF88, the 
other panels show the results adding a LER to that rate (see Table 
1). For comparison purposes, in all the panels we have drawn the 
CF88 curves as thin lines. In the panels belonging to the LER- 
1.5-10.-1.0 and LER- 1.3-0. 1-1.0 rates, the track belongs to the 
evolution of the shell that first crosses the dynamical curve. We 
assumed a chemical composition made of 50% '^C and 50% ""O 
for both curves, while the hydrostatic p-T tracks were computed 
with the realistic chemical composition obtained from full stellar 
evolutionary calculations (FRANEC code). 



drop in the carbon mass fraction at the burning place. The main 
impact of a LER on the chemical profile is a lower abundance 
of carbon in the convective core, owing to the smaller size of the 
convective region than for model CF88. On the other hand, mod- 
els igniting carbon off-center maintain the original '^C abun- 
dance at the center, because there is neither mixing nor burning 
there. However, the differences in the carbon mass fractions and 
profiles that we have found do not have a significant impact on 
the SNIa outcome. 

There is currently no successful theoretical model able to 
predict the properties of SNIa explosions univocally, starting 
from a given presupernova structure. In other words, all current 
SNIa models rely on a more or less large number of free pa- 
rameters, either related to the ignition configuration (usually in 
two and three dimensional models), or to the nuclear front prop- 
agation (in models of any dimensionality), or to both of them. 
It turns out that the explosion properties are mainly sensitive to 
these free, uncertain, parameters rather than to the details of the 
presupernova structure. In our study, the adoption of the delayed- 
detonation paradigm with a fixed transition density, pddt, and 
the use of a one-dimensional model implies that all the details 
of the presupernova are lost apart from the central density, and 
the chemical composition profile, whose most relevant byprod- 
uct is the profile of the neutron excess. As a result, the kinetic 




0^2 0.4 0.6 0,8 1 1,2 



Lagrangian Mass Coordinate (Mp) 

Fig. 4. Thermal and chemical profiles at the crossing of the dy- 
namical curves. The line types used for the models from Table 1 
are the same as those used in Fig. |2l Top panel: Temperature 
profiles. From top to bottom (at M ~ 0.2 M©) the curves 
belong to models CF88, LER-1.5-0.001-1.0, LER- 1.7-0. 1-1.0, 
LER- 1.5-0. 1-1.0, LER-1. 3-0.1-1.0, and LER-1.5-10.-1.0. Note 
that the last two models run away off-center. Bottom panel: '^C 
mass fraction profiles. For comparison purposes, we have plotted 
with a thin solid line the profile of model CF88 at the crossing 
of the ignition curve. As can be easily appreciated, the extension 
of the convective core is very sensitive to the inclusion of a LER 
in the '^Ch-'^C reaction rate. 



energy of the computed SNIa (Table O changes by less than 
~ 5% among our set of models. On the other hand, the mass of 
radioactive ^^Ni synthesized in the explosion varies by as much 
as ~ 17%. Rather surprisingly, the largest difference in the ^*Ni 
mass ejected, AM(^''Ni) ~ 0.06 Mq, belongs to two models that 
have quite similar properties during the accretion phase: LER- 
1.3-0.1-1.0 and LER-1. 5-10. -1.0. In any case, the synthesized 
masses of ^^Ni are well within the values usually accepted for 
normal SNIa, so that the presence of a LER would not change 
our basic picture of the formation of the light curve of these stan- 
dard candles. 

In Table |3] we report the yields of some elements that might 
be affected by a change in the rate of carbon burning: carbon, 
oxygen, magnesium, and silicon. We note, in particular, that oxy- 
gen is destroyed by the flame but partially regenerated in the 
outer shells of the supernova by explosive carbon burning. The 
most interesting characteristic is the monotonic dependence of 
the yield of carbon on resonance strength and energy. Although, 
in general, the variation in the carbon yield with respect to model 
CF88 is quite modest, it is largest in the two models ignit- 
ing off'-center (LER- 1 .3-0. 1 - 1 .0 and LER- 1.5-10.-1.0). However, 
the reason for the smaller carbon yield of these two models 
does not rely on the off-center location of the igniting spot. We 
have checked it by computing a mixed model in which the pre- 
supernova structure was that belonging to model LER-1. 5-10.- 
1.0, but the supernova (and the nucleosynthesis) was computed 
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Table 3. Results of the hydrodynamic explosion models plus post-processing nucleosynthesis 



Model 




M(^^Ni) 


M(C) 




M(0) 


M(Mg) 


M(Si) 






(Mo) 


(Mo) 




(Mo) 


(Mo) 


(Mo) 


CF88 


0.912 


0.379 


6.2 X 10" 


-3 


0.242 


0.045 


0.321 


LER-1.5-0.1-1.0 


0.959 


0.371 


5.3 X 10- 


-3 


0.201 


0.045 


0.341 


LER-1.5-0.001-1.0 


0.933 


0.374 


5.7 X 10- 


-3 


0.229 


0.049 


0.332 


LER-1.5-10.-1.0 


0.957 


0.351 


1.8 X 10 


-3 


0.202 


0.046 


0.357 


LER-1.3-0.1-1.0 


0.949 


0.411 


1.6 X 10 


-3 


0.189 


0.042 


0.329 


LER-1.7-0.1-1.0 


0.943 


0.369 


5.7 X 10- 


-3 


0.219 


0.045 


0.344 



Notes. Each model was computed as a delayed detonation with transition density pddt = 1.6 x 10' g cm Kinetic energy in lO""' erg. 



with the CF88 rate, the result being that the yield of carbon 
amounted to 5.4 x 10""' M©, i.e. close to the value obtained for 
the CF88 model. Hence, the reason for the small yield of car- 
bon obtained in the two models with the strongest and with the 
lowest energy resonance has to be found in the higher '^C-h'^C 
rate in the outer layers of the star, where carbon burning is in- 
complete. A higher carbon fusion rate implies a smaller flame 
width (Timmes & Wooslev 1992), which is crucial for complet- 
ing the carbon burning proces s when the flame width becom es 
comparable to the WD radius (' Dommguez & Khokhlovll201 ll) . 

Figure |5] shows the chemical profiles in the outer layers of 
models CF88 and LER-1. 5-10. -1.0 (note that the mass coor- 
dinate in this figure starts at the surface of the ejecta and in- 
creases radially inwards, i.e. to the left of the plot). While the 
chemical profile is quite similar in both models within the in- 
ner ~ 0.03 M0, carbon is burned more efficiently in the outer- 
most ~ 0.01 Mq of model LER-1. 5-10. -1.0, producing as well 
a higher yield of magnesium, neon, and silicon. Unfortunately, 
the yields of the elements most sensitive to the '^Ch-'^C reac- 
tion rate are small and are located in the lowest density regions 
of the ejecta, which implies that they can only be detected at 
a significant level during the few days following the supernova 
explosion. 

4. Physical conditions at tliermal runaway 

The structure of massive non-rotating white dwarfs at carbon 
runaway is determined by a few basic properties of the star: 
1) its central density, pc; 2) the electron mole number, Y^, or 
the neutron excess, 77 = \ - 2Y^; and 3) the number and dis- 
tribution of hot spots at which the temperature first runs away. 
We have discussed in the previous section the dependence of 
the first factor, the central density at runaway, on the presence 
of a LER. The second factor, the electron mole number, is de- 
termined by the weak interactions taking place during carbon 
simmering and by convective mixing. Even if we use a time- 
dependent mixing scheme coupled to the nuclear burning, the 
adopted nuclear network is not complete enough to describe 
the evolution of with the required precision. In addition, as 
we now demonstrate, the final value of Y^ depends critically 
on a/p. The third factor, namely the number and distribution 
of runaway centers, also known as hot spots or igniting bub- 
bles, is intrinsically multidimensional, so that it cannot be de- 
termined with our one-dimensional hydrostatic code. Hence, we 
consider alternative ways of addressing both the dependence 
of the neutronization d uring carbon simmering dChamulak et alJ 
l2008tlPiro & Bildsteni r2008) and the number and distribution of 
hot spots, on the pre sence of a LER in the '^C-n'^C reaction 
(IWooslev et al.ll2004 . 

We note that there is another potential effect of the presuper- 
nova structure on the SNIa outcome that we have not accounted 



for: the possible dependence of the deflagration-detonation tran- 
sition density on the WD chemical profile. Although this depen- 
dence has alre ady been addressed elsewhere (Bravo et al. 201^ 
IJackson et alj |2010). it is somewhat speculative and we prefer 
not to mix it with the already hypothetical character of our ghost 
resonance. 

4.1. Neutronization during carbon simmering 

The progenitor metallicity sets a ceiling on the fraction of the 
incinerated mass that can be in the form of radioactive ^^Ni, and 
thus influences the supernova lu minosity, with profound impli- 
cations for supernova cosmologv. lTimmes et al.l (i2003,) proposed 
a linear relationship between th e mass of ^^Ni and met allicity: 
M(5<'Ni) oc 1 -O.O57Z/Z0. Later. IChamulak et al.l (l2008h argued 
that even a zero metallicity progenitor is subject to neutroniza- 
tion through electron captures during the time that elapses be- 
tween the crossing of the carbon ignition curve and the onset of 
the dynamical event, i.e. until the timescale becomes on the order 
of 1-10 s. The neutronization during the simmering phase would 
affect mainly SNIa from low metallicity progenitors, and the re- 
lationship between M(^^Ni) and Z would change slightly, to be- 
come MC'^^'Ni) cc 0.965 - 0.057Z/Zo. Existing methods to de- 
termine the progenitor metallicity from the supernova properties 
depend on the degree of neutroni zation du ring carbon s immer- 
ing (,Lentz et al. 2000; Taubenberger et al.li2008; Badenes et alJ 
120081k 

The final neutron excess at runaway is a result of several 
chains of reactions taking place during carbon simmering, all of 
them primed by the '^C-h'^C reaction. The main reaction chain is 
activated by the release of a proton by the carbon fusion reaction 

'^C+'^C 2-^Na + p 

with the result that there is one electron capture for each proton 
delivered in the '^C -1-'^ C reaction. This sequence of reactions 
is fast enough relative to the fusion of '^C to ensure that every 
time a proton is released it is followed by an electron capture 
onto '^1S0 

Another significant contributor to the neutronization of mat- 
ter during carbon simmering is the electron capture onto ^''Na 

e" Na -^^^ Ne + 

that is activated also by the proton channel of the '^Ch-'^C re- 
action. However, in this case the electron capture timescale is 

- The timescale of electron captures onto '^N at the temperatures 
and densities of interest is much shorter than the WD sound crossing 
timescale, hence this reaction is always active during carbon simmer- 
ing. 
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Mass coordinate (Mg) Mass coordinate (Mg) 

Fig.5. Chemical profile in the outer layers of models CF88 and LER-1. 5-10. -1.0. The mass coordinate measures the mass within 
the white dwarf surface, i.e. the surface is at the right end of the horizontal axis. The thicker unlabeled line shows the maximum 
temperature achieved in each layer during the explosion (right axis). 



longer than the convective timescale implying that, as the tem- 
perature rises approaching the runaway, the strong interactions 
become faster than these electron captures. A direct consequence 
is that the abundance of ^^Ne freezes out at high temperatures. 

Assuming the proton yield estimated by CF88, i.e. 0.43 pro- 
tons per reaction, gives 0.14 electron captures for every '^C nu- 
clei consumed in the above chains. At densities in excess of 
10^ g cm"^, other electron captures become energetically fea- 
sible and the ratio of electron captur es per '^C nuclei destr oyed 
increases to a typical value of ~ 0.3 (IChamulak et alJl200 8'). On 
the other hand, a proton yield as low as that found by S07 for 
the 2.14 MeV resonance can modify substantially the ratio of 
electron captures per '^C nuclei destroyed. For instance, a yield 
of 0.16 protons per reaction would give 0.055 electron captures 
for every '^C nuclei consumed in the above chains, i.e. approx- 
imately a factor of three reduction in the amount of electrons 
captured for a given reduction in the carbon mass fraction. 

A reduction in the amount of electron captures activated by 
the proton channel allows other reactions to contribute signifi- 
cantly to the neutronization. The initial presence of ^^Ne opens 
up another route for electron captures activated by the alpha 
channel of the carbon fusion reaction 

'^C-h'^C ^''Ne + o' 
a-H^^Ne ^^Mg-nn 
e^H-^^Mg 2^Na, 

which has to be taken into account in the calculations of neu- 
tronization during carbon simmering. 

4.1 .1 . Impact of a low-energy resonance 

To evaluate the impact of a ghost resonance on the degree of neu- 
tronization during carbon simmering we post-processed, with 
our nucleosynthetic code, the central p-T tracks obtained with 
the hydrostatic code (Fig. [3) until they crossed the dynamic 



curve. We did not account for convective mixing, because the 
changes in the chemical composition we consider are so moder- 
ate that the contribution of fresh fuel would not have a significant 
effect on the nuclear network. 

Table |4] shows the results obtained for our reference model, 
CF88, and two different resonance strengths with Er = 1.5 MeV, 
(cjy)^ = 0.1 and 0.001 meV (first three rows). We also com- 
puted models that differ from LER-1. 5-0. 1-1.0 in terms of the a 
to proton yield ratios (last two rows). The second column gives 
dYfi/dYi2, which is actually the average, during carbon simmer- 
ing, of the ratio of the reduction in the electron mole number to 
the reduction in the '^C molar fraction. The presence of the ghost 
resonance, clearly does not affect the average number of electron 
captures per carbon nuclei destroyed unless the a or proton chan- 
nels contribute to the carbon fusion rate in significantly different 
proportions. An increase (decrease) in a/p by a factor five trans- 
lates into a ~ 50% decrease (increase) of dFe/d^n. The third 
column of Table|4]gives the average derivative of the mean mo- 
lar fraction with respect to the carbon molar fraction, which re- 
flects the progress of the nuclear network to increasingly bound 
and heavy nuclei. This quantity depends strongly on neither the 
presence of the ghost resonance nor a/p. The small decrement 
in dY\,/dYi2 for Q'/p=0.2 is due to the destruction of a third '^C 
nuclei, for each carbon fusion reaction, by a radiative proton cap- 
ture, which depends directly on the strength of the proton chan- 
nel in the '^C -h'^ C reaction. 

The main conclusion from Table|4]is the definite dependence 
of the final neutron excess on the output channel of the car- 
bon fusion reaction. In Fig. |6] we show details of the chemical 
evolution of the reference model, CF88, and the models with 
Er = 1.5 MeV and (ar/)^ - Q.I meV and different a/p. In this 
figure, the different roles played by the electron captures on '^N 
and ^^Na with respect to neutronization can be appreciated. 

The abundance of ^^Ne directly reflects the electron cap- 
tures on ^-^Na, it increases steadily at low temperatures but 
freezes out as soon as the nuclear timescale becomes smaller 
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Table 4. Neutronization during carbon simmering for different '^C+'^C rates. 



Model 


dYJdYn" 




- (Ayi2)core' 


-(AKe)co.-e'' 








(mol/g) 


(mol/g) 


CF88 


0.2182 


0.4430 


8.15 X 10^* 


1.78 X 10-* 


LER-1.5-0.1-1.0 


0.2167 


0.4170 


2.51 X lO-'* 


0.54 X 10"* 


LER-1. 5-0.001-1.0 


0.2119 


0.4232 


5.59 X 10-* 


1.19 X 10-* 


LER-1. 5-0. 1-5.0 


0.0945 


0.4517 


2.51 X 10-* 


0.24 X 10-* 


LER-1. 5-0. 1-0.2 


0.3197 


0.3920 


2.51 X 10-" 


0.80 X 10-* 



Notes. Average derivative of the electron mole number with respect to the carbon molar fraction. Average derivative of the mean molar 
fraction with respect to the carbon molar fraction. Total decrement in the '^C molar fraction within the convective core. Total decrement in 
the electron mole number within the convective core. 



than its electron capture timescale. This condition is achieved at 
r ~ 5.5 X lO** K with the CF88 rate, and at T ~ 3.5 x 10^ K with 
the (low-energy) resonant rate. After peaking at these tempera- 
tures, the abundance of ^^Ne slightly decreases as a consequence 
of the reaction 



p Ne 



,23 



Na + n . 



This last reaction, which leads to a stronger decrease in the ^^Ne 
abundance in the models as a/p decreases, does not change the 
electron mole number Hence, the peak of the abundance of ^^Ne 
gives the measure of the contribution of the electron captures on 
^^Na to the global neutronization of matter 

The evolution of the electron captures on '^N can be moni- 
tored by adding the abundances of '"*C and '^O. Following each 
electron capture on ^^N, to give '•'C, this last nucleus captures 
any available a particle and unties a chain of reactions 
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o 
o 



O-t-n 
"O + y 



Owing to the large abundance of '^O, the last neutron capture 
is only significant after the mass fraction of '^O has grown to 
X('^O) ~ 4 X IQ-'^-lQ--'. As can be seen by comparing the peak 
abundance of ^^Ne with the sum of the abundances of '*C and 
'^O, the largest contributors to neutronization in all the models 
are the electron captures onto ^^Na. 

Electron captures onto ^""Mg only contribute significantly to 
neutronizationin model LER-1. 5-0. 1-5.0, with a/p = 5. Even in 
a model with such a reduced proton yield, their contribution to 
neutronization, as measured by the peak abundance of ^^Na, is 
less than half the contribution of the electron captures onto ^^Na. 

Figure|7]shows the evolution of the decrement in the electron 
mole number as a function of the decrement in the carbon mo- 
lar fraction, for the same models as above. Each curve displays 
two parts with different slopes. The knee of the curves belongs 
to the time at which electron captures on ^■'Na freeze-out and, 
thereafter matter is less efficiently neutronized. The plot clearly 
shows that the models with the same a/p (CF88 and LER-1.5- 
0.1-1.0) are characterized by the same rate of neutronization rel- 
ative to the rate of consumption of '^C, while a different weight 
of proton and a channels implies different slopes in both parts of 
the corresponding curves. Another interesting feature in Fig. Q 
is that the total consumption of '^C in the resonant models is far 
less than in the CF88 model, and the same applies to the total 
decrease in the electron mole number. 

To evaluate the consequences of any dependence of -AFg on 
the resonance properties, we have to take into account the to- 
tal mass of '^C destroyed before reaching the dynamic curve, 
AM('^C), and the size of the convective core, Mcom (Table [TJ. 



0.7 




1 1.5 2 

-AY,, (10-^ mol/g) 

Fig. 7. Evolution of the decrement in the electron mole number 
vs. the decrement in the carbon molar fraction for models: CF88, 
LER- 1 .5-0. 1 - 1 .0, LER- 1 .5-0. 1 -5.0, and LER- 1.5-0.1-0.2. 

The total decrement in the '^C molar fraction within the convec- 
tive core, - (Ayi2)coie (third column in TablelDi is given by 



(AFi2)c 



AM(12C) 

1 2 A/rnnv 



(7) 



The total decrement in the electron mole number within the WD 
core during carbon simmering is 



(Aye)c 



d^e 

dFp 



■ (Ayi2)c 



(8) 



which is given in the last column of Table H] Accounting for 
all these factors results in a significant reduction in the neu- 
tronization when there is a LER in the carbon fusion reaction. 
The largest reduction in (Aye)coie belongs to the model assum- 
ing c/p - 5.0, in which case the reduction in the electron mole 
number in the convective core is only ~ 13% of the reduction 
obtained with the CF88 model. 

We note that the reduction in the electron mole number 
within the convective core is the result of different factors, some 
of which depend on a/p, while others depend on the resonance 
energy and strength. However, for a/p = 5, even a small strength 
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Temperature (10^ K) Temperature (10^ K) 

Fig. 6. Evolution of the chemical composition at the center of the white dwarf as a function of temperature for different '^C+'^C 
rates: CF88 rate (top left), LER-1.5-0.1-1.0 model (top right), and a/p: LER-1.5-0.1-5.0 with a/p=5.0 (bottom left), and LER-1.5- 
0.1-0.2 with Q'/p=0.2 (bottom right). To improve the readability, the upper limit to the vertical axis has been set at a mass fraction 
of 0.01, hence the curves of '^c and '*'0 (both with mass fractions > 0.1) lie out of the plot (the relative variations in their mass 
fractions are very tiny within the temperature range shown). The thick curve gives the electron molar fraction. The abundance of 
^^Ne reflects the amount of electron captures on ^■'Na: it grows quickly witii increasing temperature at low temperatures but levels 
out at r ~ 5.5 X 10^ K for the CF88 model and T ~ 3.5 x 10^ K for the rest of the models shown. The sum of the abundances of 
'^O and '"^C reflect the electron captures on '-^N: the sum grows steadily but is always smaller than the abundance of ^^Ne, except 
at very high temperatures. The top axes of the plots belonging to models CF88 and LER-1.5-0.1-1.0 are labeled with the time to 
runaway. 



of LER (for which AM(^^C) and Mconv were similar to the CF88 
model) would be enough to reduce the neufronization during car- 
bon simmering by a factor < 0.5. 



4.1 .2. Variations in either central density or accretion rate 

In reaUstic evolutionary models, the central density at carbon ig- 
nition is determined mainly by the accretion rate. However, a 
variation in M also afl'ects the thermal profile along the whole 
WD, hence the mass coordinate where carbon runs away (e.g. 
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a/p = 5.0 



12 3 4 

p (109 g cm-3) 

Fig. 8. Average derivative of the electron molar fraction with 
respect to the '^C molar fraction during carbon simmering, as a 
function of density for different carbon fusion rates: CF88 (black 
asterisks), LER- 1.5-0. 1-1.0 (blue triangles), LER-1.5-0.001-1.0 
and LER- 1.5- 10.- 1.0 (green stars), and LER-1.5-0.1-5.0 and 
LER- 1.5-0. 1-0. 2 (red and magenta circles, respectively). All the 
models shown in this panel were computed at constant density. 
The models with the extreme resonance strengths were com- 
puted only at p = 2 X 10^ 



gem 



-3 



iPiersanti et al.ll20d3h . While it is impractical to explore with our 
hydrostatic code the whole space of accretion rates together with 
all the resonance parameters, we can gain some indirect insight 
into the consequences of different accretion rates by comput- 
ing the evolution during the simmering phase with our post- 
processing code at different constant densities. Hence, we re- 
peated the calculations described in Section |4. 1 . 1 1 following, at 
constant density, the thermal and chemical evolution between the 
temperatures of the ignition and dynamic curves. Our results are 
summarized in Fig. [8] where we use four different densities in 
the range 10^- 4 x lO'' g cm"^^. 

The results show that the impact of the resonance strength (or 
even of the existence of a LER) is minimal. The points belonging 
to models with Er = 1.5 MeV and resonance strengths differing 
by four orders of magnitude are indeed indistinguishable (blue 
triangles and green stars at p = 2 x 10^ g cm"^). In contrast, the 
most influencing parameter is the ratio of the yields of protons 
to alphas in the '^C-h'^C reaction, whose impact on the variation 
in the electron molar fraction is by far higher than that of the 
density. It is remarkable that the impact of the a/p ratio on the 
electron captures is about the same at high densities as at p = 
10^ g cm even though in this case the only electron captures 
allowed are those onto '^N. 

4.2. Runaway bubbles 

The thermal structure of the white dwarf, and thus the distri- 
bution of hot spots, at thermal runaway is set during the last 
seconds before the crossing of the dynamical curve, when the 



nuclear timescale becomes of t he same order as th e convec- 
tive turnover time, ~ 10 - 100 s (IWooslev et al.ll2004 ^. hereafter 
WWK04). At this point, the timescale is so short that the tem- 
perature of the central region of the WD is no longer uniform. 
Hot spots originate, instead at random locations near the center 
and are moved by the convective cells across the turbulent core. 
During their trip, these bubbles are subject to thermal stresses 
of different sign: on the one hand, nuclear reactions release heat 
efficiently while, on the other hand, thermal exchange with the 
background cools them. If the thermal fluctuations induced by 
convection are large enough, the heat release by nuclear reac- 
tions dominates over adiabatic cooling and a bubble can run 
away in a time shorter than the turnover timescale. Although 
this process is intrinsically three-dimensional, we can model it 
following the ideas of WWK04. 

To determine the runaway temperature, T^^, of the convec- 
tive bubbles, we use the same nucleosynthetic code as before. 
Now, however, the evolution of the bubble temperature is inte- 
grated together with the nuclear network, and takes into account 
both the nuclear reactions and the convective cooling 



dT 
dT 



(9) 



where Cp is again the heat capacity at constant pressure, finuc is 
the rate of energy release by nuclear reactions including elec- 
tron captures, Sy is the rate of neutrino energy loss (although it 
is not relevant near runaway, we have included it for complete- 
ness), Vconv is the convective velocity, which we take here as a 
free parameter, and (dT/dr)^^ is the adiabatic thermal gradient. 
The last term in Eq. |9] models the loss of heat from the bubble 
to the background caused by convection (WWK04). Equation|9] 
can be easily integrated for a given initial temperature. To, and 
bubble distance to the center, ro, which we take as the radius of 
the zone that first crosses the dynamical curve. As we adopt a 
constant Vconv, we can easily calculate the position of the bubble 
at any time, f, as r = ro + Vconv?- Using these approximations, the 
initial bubble temperature. To - T^^, above which the integra- 
tion of Eq.|9]diverges and the radius, r^^, at the time of runaway 
can be computed for a given central density, po, and convective 
velocity. 

Table |5] allows a comparison between the runaway temper- 
ature, Trw, obtained with a LER and the one achieved with the 
CF88 rate. The resonance leads to a lower ignition temperature, 
from ~ 8 X 10'^ K to Trw ~ 4 - 6 x 10'^ K. The runaway 
radius, r^^, is not as sensitive to the presence of the ghost reso- 
nance, as it stays in a narrow range of r^w ^ 180 - 240 km for 
most models. The runaway radius for model LER- 1.3-0. 1-1.0 is 
larger because of the off-center crossing of the dynamical curve. 

We have explored the sensitivity of the bubble runaway 
conditions of model LER- 1.5-0. 1-1.0 to the assumed physico- 
chemical parameters: convective velocity, Vconv. density, and '^C 
mass fraction, and present the results in Table |6] As can be seen, 
for the chosen combination of parameters, the runaway temper- 
ature remains much lower than that of the CF88 model, and is 
mostly sensitive to the central density of the WD (we recall that 
the mass accretion rates that differ from those we assumed would 
lead to different central densities at runaway). On the other hand, 
a decrease in the convective velocity from 80 km s"' to 20 km s"' 
changes by only 0.2 x 10*^ K. The value of the central den- 
sity of the WD also affects substantially the runaway radius, 
which changes by a factor of two between our extreme densi- 
ties. However, it is the convective velocity that has the strongest 
influence on r^^. If Vconv is as small as 20 km s"\ the bubble runs 
away just at 90 km from the center of the WD. 
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Table 5. Runaway conditions for different '^C+'^C rates. 



Model 


P" 


X r LI 


; 

' rw 


b 


t 

'rw 




(g-cm"3) 




(10** K) 


(km) 


(s) 


CF88 


2.09 X lO'' 


0.224 


8.14 


220 


2.8 


LER-1.5-0.1-1.0 


2.33 X 10' 


0.233 


4.76 


177 


2.2 


LER-1.5-0.001-1.0 


2.38 X 10' 


0.227 


6.03 


225 


2.8 


LER-1.5-10.-1.0 


1.90 X 10' 


0.271 


4.26 


242 


1.2 


LER-1.3-0.1-1.0 


1.34 X 10' 


0.301 


4.43 


338 


1.3 


LER-1.7-0.1-1.0 


2.62 X 10' 


0.241 


5.65 


184 


2.3 



Notes. In these calculations we assumed a constant convective velocity: Vco„v = 80 km/s. Density and '^C mass fraction at the point that crosses 
first the dynamical curve. Resulting conditions at runaway: minimum initial temperature of the bubble that runs away, T^, runaway radius, r^„, 
and runaway time, ?,„. 



Table 6. Sensitivity of the runaway conditions of model LER- 
1.5-0.1-1 .0 with respect to physico-chemical parameters. 



10 



^conv 


Pc 






'"rw 




(km/s) 


(g-cm-3) 




(10** K) 


(km) 


(s) 


20 


2.33 X 10' 


0.233 


4.55 


91 


4.6 


80 


1.00 X 10' 


0.233 


5.51 


264 


3.3 


80 


3.00 X 10' 


0.233 


4.53 


157 


2.0 


80 


4.00 X 10' 


0.233 


4.26 


136 


1.7 


80 


2.33 X 10' 


0.100 


5.18 


191 


2.4 


80 


2.33 X 10' 


0.400 


4.53 


171 


2.1 



Tables |5] and |6] give the time, it takes for the bubble to 
run away after the temperature reaches T - T„. In all cases, 
the runaway time is short enough to allow for the survival of the 
hot bubbles against disruption by the Rayleigh-Taylor instability, 
which is on the order of fRx ~ 5 s (lapichino et al. 2006). The 
survival time of the bubbles against disruption sets a limit to the 
maximum distance a bubble can travel from its birthplace before 
running away. For instance, for our adopted convective velocity 



of Vrnnv - 80 km s , this distance is r„ 



400 km. 



Some examples of the temperature evolution of the bubbles 
are shown in Fig. |9] for different initial temperatures, Tq, and 
carbon fusion rates. As can be seen, the radius at which a bubble 
runs away depends sensitively on its initial temperature (i.e. on 
the level of thermal fluctuations in the convective core). As might 
be expected, the larger the excess of the initial temperature of a 
bubble over T™, flie closer it burns towards its initial location. 



4.2.1. Multispot runaway 

Once a bubble runs away, the time available for other sparks to 
ignite, 6t, is limited by the time it takes for a conductive flame 
to cross the whole region where hot spots dwell, d ~ 2rn,ax- 
As the conductive flame speed at the densi ties of interest is 
Vflame - 40 km s ' jTimmes & Wooslevll 19921) . the formation of 
multiple flame centers is possible during a time ~ 20 s. Thus, a 
multispot runaway scenario is a reasonable assumption for SNIa, 
and its properties are an important ingredient for multidimen- 
sional SNIa models. 

Because the induction of local peaks of temperature by con- 
vection is a stochastic process, we have to rely on a statistical 
approach to the multispot scenario. Indeed, WWK04 studied the 
distribution of thermal fluctuations in the convective core of an 
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Fig. 9. Evolution of the temperature of a burning convective 
bubble as a function of the distance to the center of the white 
dwarf, for diff'erent initial temperatures. To, and prescriptions 
for the '^C-h'^C rate. From top to bottom: CF88 (short-long 
dashed line), LER-1.5-0.001-1.0 (dot-long dashed line), LER- 
1.5-0.1-1.0 (solid line), and LER-1.5-10.-1.0 (long-dashed line). 
For each reaction rate, we followed the evolution of the bubbles 
starting at the center, r = 0, with different trial initial temper- 
atures, until we found the minimum initial temperature leading 
to runaway of the bubble before reaching the end of the con- 
vective core. The distance to the center varies linearly with the 
elapsed time, r - + Vconvf, where Vconv = 80 km s"' is the 
convective velocity. To test the impact of the off'-center ignition 
(at /;,w ~ 149 km) in model LER-1. 5-10.-1.0, we computed the 
bubble evolution with the same rate but starting from the center. 
The runaway temperature, shown by the thin long-dashed line, 
is nearly the same as in the off'-center calculation of the same 
model. 



igniting WD and proposed that there is an exponential probabil- 
ity density function (EPDF) of temperatures 



^ cc-(l - f)ff-r»)/Ar 

dr J' 



(10) 
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where / is a parameter correlated to convective mixing effi- 
ciency, which following WWK04 we assume to be / ^ 0.9, 

is the average temperature, which we take without loss of 
generality as the runaway temperature, T^^ - T„, and AT ~ 
0.1 X 10** K or sHghtly lower (WWK04). To derive a distribu- 
tion function for the runaway radii (distance to the center of the 
WD) of the bubbles, we need the relationship between and 
the initial temperature excess over T^^- We computed, for both 
model CF88 and model LER- 1.5-0. 1-1.0, the temperature evo- 
lution of bubbles with an initial temperature, Tq, above the run- 
away temperature, from which we found the coiTesponding run- 
away radii, and we fit an empirical relationship between To and 

(Fig. [Tol l. A second-order polynomial provides a good de- 
scription of (To) and allows us to find the probability density 
function of runaway radii 



dP 
dr,-w 



dP/dT 
dr„/dTo 



(11) 



The results are shown in the bottom panel of Fig. [10] The distri- 
butions of runaway radii for both carbon fusion rates are quali- 
tatively similar in shape and width, but with the LER the peak of 
the distribution is ~ 50 km closer to the center of the WD. Such 
a displacement of the runaway radius can have important con- 
sequences for SNla models such as the gravit ationally confined 
detonation model (see, e.g.. lRopke et al.ll2007h . 



5. Conclusions 

We have systematically explored the consequences that a hypo- 
thetical low-energy resonance in the carbon fusion rate would 
have for the physics and outcome of thermonuclear supernovae. 
We have considered resonance energies in the range ZiR = 1 .3 - 
1.7 MeV, with strengths limited by the available experimental 
cross-section data at E^m > 2.10 MeV. We have also studied 
different ratios of the yields of a to protons, a/p - 0.2 - 5.0. 
Within these limits, the phase of the supernova that is most af- 
fected by the presence of a LER is carbon simmering. In partic- 
ular, the degree of neutronization of the white dwarf prior to the 
explosion can be reduced by a large factor even for a LER with 
a strength much smaller, say by four-six orders of magnitude, 
than that measured for the resonance with the lowest experimen- 
tal energy (Er = 2.14 MeV, S07). 

The main changes induced by a LER with respect to the 
model that uses the standard carbon fusion rate, CF88, are: 

- The central density of the white dwarf at runaway can 
increase by ~25-30%, for a given mass accretion rate. 
However, a strong enough resonance gives an off-center ig- 
nition, in which case the density is lower than that obtained 
with the CF88 rate. 

- The minimum neutron excess of thermonuclear supernova 
ejecta, determined by the progenitor initial metallicity and 
by several electron capture reactions synchronized with car- 
bon simmering, is quite sensitive to the relative a to proton 
yield of the carbon fusion reaction The presence of a small 
resonance at 1.5 MeV (1.3 MeV) with a strength > 13 neV 
(> 0.12 neV) and an alpha/p yield ratio substantially differ- 
ent from unity would dominate over the non-resonant con- 
tribution at the relevant Gamow energies and, hence, would 
have a non-negligible impact on the neutronization of matter 
during carbon simmering. 

- The mass of the convective core is also affected by a LER. In 
general, the lower the resonance energy and/or the stronger 
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Fig. 10. Distribution of runaway radii. Top. Runaway radius as a 
function of the temperature above the runaway temperature (T™, 
see Table[5]): stars belong to model CF88, and triangles to model 
LER- 1.5-0. 1-1.0. The lines are quadratic fits to the dependence 
of runaway radii on temperature: r - ax^ + bx + c, where r is in 
km and x = ln (Tq - T^), with both temperatures in 10^ K. The 
parameters of the polinomial are a = -2.87, b = -64.08, and 
c = -13.09 for the CF88 rate, and a = -2.72, b = -60.59, and 
c = -45.21 for the LER- 1.5-0. 1-1.0 rate. Bottom. ProbabiHty 
distribution function of runaway radii, in arbitrary units. The 
probability distribution function has been obtained as the prod- 
uct of the exponential probability distribution function (EPDF) 
of temperatures proposed by WWK04 (eq. 38), dP/dT, and the 
inverse of the derivative dr/dT computed from the polinomial 
fit to r(T). The curves belong to the CF88 rate (short-long dash) 
and the LER-1.5-0. 1-1.0 rate (solid). 



the resonance, the smaller the convective core. The neutron- 
ization produced during carbon simmering erases the im- 
print of the initial metallicity of the supernova progenitor, 
and convection homogenizes the core composition before the 
thermal runaway. Hence, the smaller the convective core the 
higher the mass of the supernova that bears a record of the 
initial metallicity. 

- The time it takes the WD to cover the path from the ignition 
curve to the dynamic curve (simmering time) is substantially 
shorter (by up to a factor ten) if there is a strong LER. 

- The runaway temperature is strongly affected by the pres- 
ence of a LER, which could lower by as much as 2 - 
4 X 10** K. 

- In a multispot runaway scenario, the distribution of ignit- 
ing bubbles changes quantitatively in the presence of a LER. 
With such a resonance, even of moderate strength, bubbles 
run away closer t o the center than predicted by using the 
CF88 rate (see also llapichino & Lesaffrell20Tol) . On the other 
hand, if the resonance were strong enough to provoke an off- 
center ignition, it would have the opposite effect. 

- At the present level of theoretical understanding of the 
physics of thermonuclear supernovae, it is difficult to pre- 
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diet the true impact of an increase in the rate of carbon fu- 
sion on the nucleosynthesis and energetics of the explosive 
event. Such an study would require a completely parameter- 
free successful SNla model, which is currently not available. 
Within the one-dimensional delayed-detonation paradigm, 
we have found quite small variations in the kinetic energy 
of the supernova, K, mild variations in the yield of ^^Ni (but 
still well within the observationally allowed range), and a ra- 
tio of K to the mass of unburned carbon ejected that is quite 
sensitive to the existence of a LER. 

Finally, we briefly discuss some potential (even though spec- 
ulative) implications of the distributions of the runaway radii of 
the ignition bubbles that we have found for the different mod- 
els with/without a LER. We have shown that both the runaway 
radius, r„, is affected by a LER (Table |5]l, and the runaway 
bubbles cluster around (Fig. [TOl i. However, although a LER 
may influence the outcome of the explosion, it is impossible at 
present to make a more specific prediction. On the other hand, 
we have shown that r,„ is a sensitive function of the central den- 
sity (Table|6]), with the result that at a higher density the igniting 
bubbles would be concentrated toward smaller runaway radii, 
which can have important consequences for the subsequent de- 
velopment of the explosion. The effects of the distribution of 
hot spots in the WD at runaway has been explored just recently 
with the aid of multidimensional SNla models. For instance, 
ISeitenzahl et alj ( 1201 II) searched for the systematic effects of the 
central density on SNla luminosity by simulating WD explo- 
sions starting from hot bubbles randomly located within a radius 
independent of p^. As we have shown, it would be advisable to 
use a density-dependent distribution of the bubble positions in 
these studies, even if the carbon fusion rate is that of CF88. 

We conclude that a robust understanding of the links between 
SNla properties and their progenitors will not be attained until 
the '^C-(-'^C reaction rate is measured at energies ~ 1.5 MeV. 
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